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Abstract 

The multiparton density matrices of the QCD gluon cascade are investigated. 
The generating functional and master equation for momentum space multiparton 
density in the double logarithmic approximation (DLA) are proposed. 

1 Introduction 

Recently, several groups have analyzed in great detail the multiparton distri- 
butions in the QCD gluon cascades The results of the investigations show 
that perturbative QCD [0 provides a powerful framework not only for the 
description of hard quark and gluon jets but also of much softer multiparti- 
cle phenomena. Although not understood, the hypothesis of parton-hadron 
duality provides an apparently successful link between theoretical parton 
distributions and observed particle spectra. This prescription was exten- 
sively tested in single particle spectra ( and total multiplicities) and found 
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in a good agreement with available data (see e. g. 0]). Recently, there ap- 
peared indications that it may also work for multiparticle correlations [^]. 
These unquestionable successes invite one to study further consequences of 
the theory for processes of particle production. At this point we would like to 
notice that, if one wants to exploit fully the quantum-mechanical aspects of 
the QCD cascade, it is necessary to study the multiparton density matrix. 
The multiparticle distributions calculated so far give only diagonal terms 
of the density matrix and thus represent rather restricted ( although very 
important ) part of the information available from the theory. 

It is perhaps important to stress that in contrast to what is usually be- 
lieved, the interest in studying the multiparticle density matrix is not purely 
academic. As suggested in density matrix allows to obtain the multipar- 
ticle Wigner functions and consequently gives information about the space- 
time structure of the system. This allows to make predictions on the shape 
and range of the HBT interference with clear experimental consequences. 

We present in this paper a study of the multiparton density matrices of 
the QCD gluon cascade in the double logarithmic approximation (DLA) [|| 
of QCD. We propose a master equation for the functional which generates 
momentum space density matrix of arbitrary order. The paper is organized 
as follows. In Section 2 the definition of the density matrix is given. In 
Section 3 we recall briefly the main assumptions of the DLA approximation 
for the QCD cascade, and the generating functional method for multiparticle 
densities. In section 4 we propose a scheme to calculate the density matrix in 
DLA formalism. We formulate the master equation for the functional, which 
can generate density matrices at first at the most general level, and later in 
the approximation of the "quasi - diagonal" limit. The approximation allows 
to derive explicit integral equations for density matrices of arbitrary order. 
To show explicitly, how the method works, we derive master equation for 
single and double inclusive particle density matrices in the last section. 
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2 Definition of the density matrix 



We would like to investigate the properties of the momentum space density 
matrix. Let us consider any particle production process where all the pro- 
duced particles are already "real" i.e. they are on mass shell. Then the 
exclusive m-particle density matrix takes the form : 

m 

d^'^iK^ • ■ ■ , k;,; ki, . . . , k^) = S{D'){K, k'jS^D){ki, ...,km) n(4a;fc^'^ikJ~ 

D',D i=l 

(1) 

where S'(i5)(/ci, . . . , km) denotes the probability amplitude to obtain exactly 
m particles with the momenta ki, . . . ,km ior the given Feynman diagram D, 
and the sum is taken over all possible realization of the process i.e. over all 
possible Feynman diagrams for each amplitude separately. We have taken 
into account also phase space factors. 

The inclusive m-particle density matrix is then given by : 

/ . 7 rr / "^m+if^ln d^^i)^\, . . . , \s! , k^+i, . . . , k„; ki, . . . , k„), (2) 
^ (n-m)! J 

where di[k]j = (Pki . . . (fkj. 

The "diagonal" elements of the exclusive and inclusive density matrices 
give the multiparticle exclusive and inclusive densities respectively : 

p^'^(ki, . . . , k^) = d^%ki, . . . , k^, ; ki, . . . , k^) (3) 
P (ki, . . . , k^) = d (ki, . . . , k^, ; ki, . . . , k^). 

The on-mass-shell Fourier transforms of the density matrices define the 
exclusive and inclusive multiparticle densities in the configuration space. 
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3 Cascade in momentum space. DLA formal- 
ism 

In this section we recall the main assumptions of the DLA approximation in 
momentum space 0. This approximation accounts for the leading double 
logarithmic contributions to multiparton cross section, but it describes quite 
well the structure of the parton cascade at high energies. We will shortly 
review the formalism for the simplest QCD process: e~^e~ qq + g . . . g. To 
simplify the analysis one uses here the planar gauge. In this gauge q and q 
emit soft gluons independently since the interference between two emission 
amplitudes vanishes. Virtual corrections appear in the form of the QCD form 
factors. 

3.1 Multigluon amplitudes 

For any tree Feynman diagram D ( without 4-gluon vertices) the amplitude 
to obtain exactly m gluons (see Fig. 1 ) emitted by the initial q [q] of 4- 
momentum P with the account of the virtual corrections reads : 

Se....e^{K ...,k^) = i-ire-"^ n e-"^, (4) 

1=1 

and Mp^{ki) equals : 

MpXh) = gs^^-^^QpMGp^, (5) 

where: 

gs = Vena's, 

n is the number of gluons emitted of quark (antiquark), 

ki = {uji, kj) denotes the 4-momentum of the i-th soft gluon, 

Cj = e'f\ j = 0, . . . ,3 describes its polarization 4- vector, 

Pi is the 4-momentum of the parent of the ith gluon, 

Gp. denotes the color factor for the given vertex of the tree diagram D, 
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Qp.{ki) is a generalized step function which restricts the phase space for ki 



ep(fc) : = a; < P°, ^kP < 0, cjOy^p > Qo}, (6) 

where P is the momentum of the parent of a given parton k, 9 denotes the 
emmision angle of the previous parton splitting on the P— line, and Qq is a 
cut-off parameter. 

Radiation factor in (4) contains the function w{Pi), defined as: 

w{P) = / d'k I Ap(k) \l, (7) 
where Ap{k) introduces the phase space factor in the form : 

The expression | Ap{k) \l in (7) has been already averaged over two physical 
transverse polarizations e^,e^. It should be also emphasized that produced 
gluons are "real" (on- mass-shell) particles, so the energy cjk of the gluon of 
the momentum k can be approximated as : 

u;k^\k\. (9) 

Summing of the color factors G over the color indices gives the result : 

where P denotes the 4-momentum of the quark ( antiquark) which initalizes 
the gluon cascade. 

3.2 Multigluon densities 

In the DLA approximation different tree diagrams come from different non- 
overlapping kinematic regions. Hence they do not interfere. Therefore, to 
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calculate exclusive and inclusive multigluon densities it is enough to sum up 
the amplitudes (4) incoherently over all possible tree diagrams D. Hence one 
gets exclusive density pf? (ki, . . . , k^) in the form : 

m 

p^^(ki,...,kJ = ^nMp.(^^)le. (11) 

D i=l 

parametrized by the momentum P of the quark (antiquark) which initializes 
the cascade. Multigluon inclusive density pp (ki, . . . , k^) follow from the 
formula (2), (3) : 

p-(ki, . . . , k„) = ^ _ / p^^(ki, . . . , k„). (12) 

n=m v'' •' 

3.3 Generating functional 

Now the difficulty is how to perform the summation in (11) and (12) over 
all diagrams D in the convenient way. The problem has been solved by 
introducing the method of generating functional (GF) (see and references 
therein). Generating functional for multigluon densities Zp[ii\, which fulfills 
the master equation : 

Zp[u] = e-'^(^) exp{ J I Ap{k) \l u{k)Zk[u]), (13) 

reproduces contributions of all possible tree diagrams D and allows to express 
multigluon densities pp (ki, . . . , k^) and pf? (ki, . . . , k^) as : 

p^^(ki, . . . , k„) = — Zp \{u=o} (14) 

OUi . . . OUm 

and 

p^p"(ki, . . . , k„) = ^Zp (15) 

OUi . . . OUm 

where Ui denotes u{ki) and respectively ^ = j^^- 

We would like to emphasize once more here the simplicity of description 
of the multigluon densities (14) and (15) in the language of the GF. The 
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method allows to forget the complicated summation procedure especially in 
the case of inclusive densities, and express the multiplicities in a simple, 
compact form (for details see [§]). 

4 Calculating of the density matrix in the 
DLA formalism 

The DLA formalism in momentum space gives a good description of the struc- 
ture of the gluon cascade. The GF scheme suggests a clear recipe how to 
construct multigluon densities, and allows to describe and investigate prop- 
erties of the gluon distributions in a very convenient way. 

In this chapter we would like to discuss calculating of the density matrices 
in the DLA approximation for the same process namely : e^e~ —>■ qq + 
g . . . g. First we derive the general expression for the exclusive and inclusive 
density matrices d^p{k[, . . . , k^; ki, . . . , k^) and d^p(k[, . . . , k^; ki, . . . , k.^). 
The task looks more complicated because in this case the interference between 
different diagrams in (1) generally does not vanish. Let us define two GF, 
Zp[u] and Zp,[w], which generate the sum of all the tree amplitudes (4) and 
the sum of their complex conjugates respectively : 

Zp[u] = e-""^^^^^ expi J £kAp{k)uik)Zk[u]) 
Z*p,[s] = e-"'^^')/' exp( J d^kA*p,{k)s{k)Z;[s]), (16) 

The multigluon density matrices can be then expressed as (see Appendix A) : 
dpiK, • • • , k^; ki, . . . , k„) = — Zp[u]Z*p,[s] \{u=s=o},p=p' 

OSi' . . . OSm' OUi . . . OUm 

Using the formulae (17) one can obtain complicated equations for the 
density matrices. The summation problem does not dissapear. It is only 
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hidden in the compact form of (17). However, in the above formula we have 
taken into account all the possible interferences between different trees D and 
D'. Detailed analysis of the DLA gives the result that not all the diagrams 
mix up: one can distinguish classes of interfering diagrams. Nevertheless, at 
the general level we did not succeed in formulation of such a GF which would 
include only these interfering diagrams. 

However, we do not need the most general formula for the density matrix. 
We are interested in its behaviour if the differences | ki — k[ | , . . . , | km — k'^ \ 
are small ( large momentum differences will not contribute to Fourier trans- 
forms). It can be checked that in this limit interferences between different 
diagrams vanish, and one can sum up only "squared" contributions from 
the same graphs. So we construct the new GF for density matrices which 
reproduces only these relevant contributions of the tree diagrams and their 
virtual corrections. The last step will be to derive explicit integral equations 
for single and double particle inclusive density matrix. 

4.1 Density matrix in quasi- diagonal limit 

We investigate the density matrix in the so-called quasi-diagonal limit i.e. if 
for the m-particle matrix the differences | fci — /c'j^ | , . . . , | /c^ — A;^ | are small. 
This requirement simplifies the problem substantively. From the analysis of 
all diagrams contributing to single particle density matrix d'^p{k[, ki) it can 
be proved (see Appendix B) that interference between different diagrams 
appears only if | a;i — a;i/ |>> or | 9ip — ^i/p |>> 0. This statement can 
be generalized for any m-particle density matrix. If we have m particles, and 
(m-1) from them are "close" to each other e.g. ki = k[^ . . . , km-i = ^m-i 
then the interference of the different diagrams will take place only if either 
energies or angles of km and k'^ are strongly ordered. 

So, in our approximation we can exclude the interference between the 
different diagrams, and sum only the "square" contributions from the same 
ones. Hence, the exclusive and inclusive density matrices can be expressed 
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as 



D i=l 

< Sei...em.{K-> • • • ) k'jn)Se^...em{kl, • ■ ■ , km) >(ei...em)) (l^) 

oo -1 . m 

4"(ki, . . . , k:„; ki, . . . , k^) = ^ _ Yl / dm+i[k]n UiAcOkicokr'^' X 

n=m m). ^ J 

n 

X [| (4a;jfc.a;fe.)~^/^ < 5ei...e„(^'l, • • • , C> ■ ■ ■ > kn)Se,...er,{ki, ■■■,kn) >(ei...e„)(l9) 

j=m+l 

The summation in (18), (19) over D can be easily performed using the 
GF we define in the next subsection. 

4.2 Generating functional for density matrices in the 
quasi- diagonal limit 

We propose a master equation for the new generating functional Zpfp[u{k', k)] 
in the form : 

oo 1 „ 

ZMu] = e-^(^' ^) J2 -1 / d,[k%d^[klu{k[, k,) . ..u{k'^, k^) x 

X < A*p,ik'^)Ap{ki) >ei . . . < A*p,{k'JAp{kn) >e„ ^fc^fciM • • • ^fc^fcj^] X Pl',...,n';l,...,n,(20) 

where the function -Pi',...,n';i,...,n introduces the requested parallel angular 
ordering for n particles in the form : 

Pl',...,n';l,...,n = ^{^k[^P' > ■ ■ ■> ^k[^P') ^{^ki^P > 

(ii,...,i„)ePerm(l,...,n) 

and the radiation factor fulfills the condition : 
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kinP)' 



(21) 
(22) 



The multiparticle density matrices can be then expressed as (for proof see 
Appendix C) : 

(^1, • • • , k ; ki, . . . , km) = Zpip I r„=o},p=P' (23) 

bUm',m ■ ■ ■ OUv^i 

dp (k'l, . . . , k^; ki, . . . , km) = Zpip \{u=6^{i'-i)},p=p' ■ (24) 

OUm',m ■ ■ ■ OUy^i 

If the profile function u{k' , k) equals to and 5^{l' — I) respectively the 
functional Zpip takes the form : 

Zpp[u = 6%l' - 1)] = 1. (25) 

The GF from (20) is the generalization of the GF defined in (13). Using 
GF (20) one can also easily obtain multigluon densities (11) and (12) in 
momentum space : 

Pp (ki, . . . , kjn) = Zp'p |{„=o},p=P' (26) 

dUm,m ■ ■ ■ Sui^i 

in ^"^ 
Pp(ki, . . . ,km) = T ^ Zp'p \{u=SS{l'-l)},P=P' ■ 

For further analysis wc need the exact value of the product of the am- 
plitudes < A*p,{k')Ap{k) >e which occurs in (20). Averaging over gluon 
polarizations (see Appendix D) gives the result : 

< A*p,{k')Ap{k) >e= Ap,p{k',k) = AglG*p,Gp^l==-—^Qp,{k')Qp{k) 

y Auj^'^ixJ'^ ypk"p'k' 

(27) 

4.3 Single and double inclusive density matrix 

From the formula (24) one can derive iterative equations for inclusive density 
matrix of an arbitrary order. As an example we present equations for single 
and double particle density matrices : 
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and 



d']i{k[] ki) = J d^kApp{k, k)di^{k[; ki) + 
+fp{k[, k,)App{k[, k,) Zr,,[u = S%1' - I)] (28) 

d'^ik'^, k'^- ki, k2) = J d^kApp{k, k)d'^{k[, k'^- ku k2) + 

+App{k[, = 6%l' - l)]fp{k[, fci) + 

+App{k'„ k2)-^Z2',2[u = d'\l' - l)]fpik'„ k2) + 

+App{k[, h)Zy,[u = - l)]App{k'^, k2)Z2',2[u = - I)] fp{K, ki, k2) + 
+App{k[, ki)Zri[u = - /)] / d^'kAppik, k)d^^{k'^, k2)fp{k[, k; ki, k) + 



+App{k'^, k2)Z2'2[u = - I)] J d^kApp{k, A;)4"(A;;, k,)fp{k'^, k; k2, k) + 

+ {j d'^kApp{k,k)dt'{k[,k,)) {J d^kApp{k,k)dt'{k'^,k2)) (29) 
Factors fp{k[, . . . , k'^; ki, . . . , km) can be expressed as : 

fp{k[, . . . ,k'^]ki, . . . ,km) = e~'^(^'^)^— / di[p]nApp{pi,pi) . . . App{pn,Pn) x 

n=0 



X-Pl,...,n,fci,...,fc;,;l,...,n,fei,...,A:™(30) 

For m = 1 we have calculated fp{k', k) explicitly. It equals to : 

fp{k', k) = Q{du'P > efep)e^«^^(^" ^-"^ ^) + Q{dk>p < ekp){k' ^ k) 

(31) 

Equations (28), (29) in the limit k[ = ki, k2 = k2 reduce to multiparticle 
density equations ||^ as expected. However, also for any k, k' the structure 
of (28), (29) looks quite similar to the structure of multiparticle density 
equations. This suggests to get the solution using the technique similar to 
the technique described in M. 
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5 Conclusions 



Our conclusions can be listed as follows : 

1. We have obtained master equation for the density matrices GF. Prom 
this relation one can derive equations for density matrices of arbitrary order. 
Integral equations for single and double inclusive densities look quite simi- 
lar to the equations for multiparticle densities we already know. It allows 
to expect that they can be solved in a way similar to that used for multi- 
particle densities. Hence, as the next step in our approach we will try to 
derive multiparticle density explicitly and to investigate their behaviour in 
the coordinate space. 

2. We know that the GF "works" in the quasi-diagonal region but we 
do not know what it really means quantitatively. The further quantitative 
analysis is needed. If the quasi-diagonal region is large enough, we could 
easily obtain the multiparticle densities in the configuration space taking the 
on-mass-shell Fourier transform of the density matrix (it needs also to be 
checked if the functions we want to transform are "smooth" enough to get 
the Fourier integral vanishing in large \ ki — k[ \, . . . , \ — k'^ \ limit). 

3. One can try to find the general GF for the density matrices which would 
work also far from the quasi-diagonal region. It would solve the problems 
mentioned above. Nevertheless, it seems to be quite complicated. Till now, 
we did not find the rule which could iteratively distinguish classes of the 
interfering diagrams. 
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6 Appendix A 



Punctionals Zp[u] and ^p'[s] defined in (16) generate the sum of the am- 
phtudes (4) and their complex conjugates over all possible tree diagrams 
respectively. It can be seen when one rewrites i.e. the master equation for 
Zp in the form of the diagram series (see Fig. 2). Prom that construction 
follows the form of the exclusive density matrix (17) : 

dp{K, • • ■ , k;,; ki, . . . , k„) = -z J — Z*p[s] \{s=o} T 7 — Zp[u] |{„=o} 

dSi' . . . dSm' OUi . . . dUm 

(32) 

Then the inclusive density matrix calculated from (2) and (32) looks like : 



X 



^^Q'"- " 5sSu 5sii...SSr 

^— -^ZpM |,„.o} (33) 

and from the identity for the product of any two functionals F, F' : 

6 5 

n=0 



F[u]F'[w] |,..^,..o}= E ^ (/ d'kj^s^m^] ku^oy F'[w] |,_o} (34) 
follows (17). 



7 Appendix B 

In DLA one gets for the two gluon contributions four different tree graphs 
Ma, Ml,, defined on non-overlapping kinematic regions (see Fig. 3 ). 

Emmited gluons are either angular (AO) or energy ordered (EO). 

Let us consider all the diagrams contributing to the single particle density 
matrix d'^p{k[,ki). From the (AO) and (EO) follows that the interference 
between any two different graphs will appear only, if either energies uji , u[ or 
emmision angles 6ip, Oyp of produced gluons are strongly ordered (Fig. 4). 
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This statement can be generalized for any m-particle density matrix by 
induction. If we have m particles, and (m-1) from them are "close" to each 
other eg. ki = k[, . . . , km-i = k'^-i then the intereference of the different 
diagrams will take place only if either energies or angles of km and k'^ are 
strongly ordered. 



8 Appendix C 

Zp>p[u{k', k)] defined in (20) produces correct exclusive density matrices (23). 
This statement follows from the master equation for Zp'p represented in the 
form of diagram series (see Fig. 5). The series reproduces all "squared" 
contributions of the same tree diagrams, and excludes interference of the 
different ones (function -Pi',...,n';i,...,n)- 

The form of the inclusive density (24) follows from the formulae (2) and 
(23). Substituting (23) into (2) we get as a result : 

°° I r () 
d^{k[, ...,k'^;k„...,kJ = Y.-( d?kd?k'8\k' - k)- T X 

X-r -r Zp,p[u] |{„=0} (35) 

dUm>,m ■ ■ ■ OUv^i 

which represents the functional Zpip[u] expanded around "null" in the 
"point" u^5^{l' -I). 



9 Appendix D 

We want to average relation (27) over the physical polarizations of the pro- 
duced gluon. The exact expression to be summed over polarizations looks 
like : 

{k-P){k'-P') ^ ' 

where e, e' are polarization of the same gluon (one gets e' from e taking the 
limit k' = k). The gauge fixing we use in the approach allows to neglect 
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contributions to (36) coming from the " nonphysical" polarizations e° and e^. 
Furthermore, it requests time components of the physical polarizations e^, 
to be equal to 0. Space components of e^, can be then constructed in the 
form : 

P X k 



6^ = 



P X k 



and respectively for e' 



2 P' X k' 

e' = -, 

P' X k' 



= (38) 
I C'^ X k' I 

Summing over these 2 polarizations and expanding scalar products in (36) 
one obtains finally : 

(e^ -P)(e'^ -FO ^ ^ ^_ 
{k ■ P){k' ■ P') ujuj' OpkOp^k' ^ ' 

The result can be easily confirmed for any two physical polarizations e^, 
lying in the plane e^e^ : 

e"*" = cos^p + sin^p 

= —sinipe^ + cosipe^ (40) 
and respectively for e'^, e'^ with the phases (p — ip'. 
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